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Abstract. In the poster presented in Cool Star 15, we analyzed the effect of disk accretion on the 
, evolution of very low mass pre-main sequence stars and young brown dwarfs and the resulting 

■ uncertainties on the determination of masses and ages. We use the Lyon evolutionary 1-D code ((Qj]) 
assuming a magnetospheric accretion process, i.e., the material falls covering a small area of the 
radiative surface, and we take into account the internal energy added from the accreted material 

■ as a free parameter e. Even if the approach to this problem is phenomenological, our formalism 
provides important hints about characteristics of disk accretion, which are useful for improved 

O ■ stellar interior calculations. Using the accretion rates derived from observations our results show 

l-j ' that accretion does not affect considerably the position of theoretical isochrones as well as the 

luminosity compared with standard non-accreting models (e.g., |2j]). See more discussions in a 

forthcoming paper by Gallardo, Baraffe & Chabrier (2008). 
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^ : INTRODUCTION 

m ! 

! The origin of accretion process is still not well understood, but is it clear that it is strongly 

linked to the evolution and lifetime of the surrounding circumstellar accreting disk. It 
appears likely that, in the T-Tauri phase, the inner disk is disrupted by the magnetic field 
of the central object in which the disk material is channelled along magnetic field lines 
to crash onto the star ([HI]). The accreted material falls onto the star in accretion ring or 
spot, assuming that magnetospheric accretion column covers only a small fraction of the 
^ . stellar photosphere. 

One of the main physical quantity controlling accretion material through the disk is 
the mass accretion rate, M. Recent observations of young low mass stars and brown 
dwarfs formation regions derived accretion rates increasing with the stellar masses: 
M oc Ml flUl]; IH). On the other hand, the derived sub and/or stellar masses depend 
on evolutionary tracks. However, these evolutionary models at young ages suffer from 
several uncertainties, mainly due to unknown initial conditions resulting from prior 
cloud collapse and accretion phases. The evolutionary tracks need to be modified to 
take into account accretion process that build up stars, 

The purpose of this poster is to present theoretical evolutionary models analyzing the 
effect of accretion on the evolution of young very low mass pre-main sequence stars 
and young brown dwarfs and the resulting uncertainties on the determination of masses 
and ages, based on detailed input physics describing the structure and evolution of such 
low mass objects and to extend it to the brown dwarfs regime. One of the main idea is 



to analyze the sensitivity of evolutionary tracks to different accretion rates in the stellar 
sub-stellar mass regime. 

RESULTS 
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FIGURE 1. Mass-radius relationship for an accreting object with an initial mass of 0.1 solar masses, 
constant accretion rate of dM / dt=5E-7 M Q and different £ values (showed in the inset). The polytropic 
model corresponds to the [6] formalism. The internal energy U is added to the accreting object as 
AU « eAM, where M is the object's mass. 
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FIGURE 2. Hertzsprung-Russell diagram (HRD) for an accreting object with initial mass of 0.05 Mq 
and e = 0.0 (i.e., without internal energy added to the object) represented in solid line. The vertical short 
dashed lines are cooling tracks for non-accreting low-mass objects, with masses indicated near the curves 
(in Mq). The full squares indicate the position of non-accreting objects with the same age (indicated 
by the numbers close to the squares, in Myr) and same mass than accreting counterpart (indicated by 
the triangles just below the squares). The numbers in brackets (close to the triangles) give the age, in 
Myr, of non-accreting objects at the position indicated by the triangles. We see that assigning an age 
for an observed young object with a given luminosity from non accreting tracks, we can significantly 
overestimate its age. This illustrates the uncertainty in age determination based on standard evolutionary 
tracks at young ages. 
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FIGURE 3. Effect of accretion on the mass-radius relationship. Solid line shows an accreting object 
with initial mass of 0.05 M Q . The accreting model was calculated for e = 0.0. The short dashed line 
indicates the radius of non-accreting objects with the same mass and the same age as their accreting 
counterparts. Ages for the accreting object, in Myr, are indicated by the numbers close to the squares. 
At any time the structure of the accreting object is more compact than its non-accreting counterpart, i.e, 
accretion time-scale becomes of the order of Kelvin-Helmholtz time-scale, thus the accreting object does 
not have time to expand to the radius at the corresponding mass it would have in the absence of accretion. 



P- 



: .5 



0.5 




0.6 



0.1 0.2 0.3 0.4 0.5 

FIGURE 4. Mass-radius relationship for an accfetirlg object with initial mass of 0.05 M & and an 
accretion rate of 5e-7 M & /yr (for all showed models). The solid line displays a model using £ = 0.0, 
dotted long dashed line shows a model for e — 0.1, and long dashed line portrays a model with e = 0.5. 
The short dashed line correspond to the short dashed line in Fig. 3. Following the discussion by 0], the 
accreting object expand in radius and it evolves maintaining a nearly fixed mass-radius ratio as long as 
deuterium (D) burning nuclear energy is sufficiently high. The effect of e ^ makes the star grow in 
radius until the extra supply of D becomes ineffective, thus the luminosity loss is dominant and the object 
need to contract. 
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FIGURE 5. Evolution of deuterium abundance for an accreting object of 0.05 Mq. Solid line shows 
a non-accreting model. The dotted long dashed and short dashed lines are models calculated using 
dM / dt=\E-% M@/yr and e = 0.0 and e = 0.1 respectively. The long dashed and dotted lines are models 
using dM / dt=lE-6 Mq /yr and £ = 0.0 and e = 0. 1 respectively. The additional internal energy compen- 
sates the energy loss by radiation at the surface, because the internal energy input reduces the need for 
deuterium burning to supply the entire luminosity. 
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FIGURE 6. HRD for young low mass objects and several theoretical isochrones taking into account 
different accretion rates and £ values. Data were taken mostly from JUl, and J3l showed in filled 
circles. Standard non-accreting isochrone from 0.01 to 0.1 Mq (indicated by the numbers close to the 
squares) is displayed in solid line. The long dashed line shows a model using dM / <if=lE-8 M cl ,/yr and 
£ = 0.0. Dotted line shows a model using dM / dt=5E-7 M & /yr and £ = 0.3. Short dashed line corresponds 
to a model using dM / dt=lE-l M & /yr and £ = 0.5. All models were calculated for 1 Myr old. Data age 
range from 1 to 3 Myr old. 
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